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HANEX  SIMULATIONS  OF  LASER-TARGET  INTERACTIONS: 

I.  RADIATION  PROCESSES 

I.  INTRODUCTION 

Over  the  past  several  years  the  Plasma  Physics  Division  at  the  Naval 
Research  Laboratory  has  run  a  laser-target  experiment  with  the  objectives  of 
modeling  physical  processes  occuring  in  a  high  altitude  nuclear  explosion 
(HANE)  This  latter  program  is  under  the  auspices  of  the  Division  of 
Atmospheric  Effects  at  the  Defense  Nuclear  Agency.  In  this  regard  the 
experiment  is  intended  to  provide  a  physical  understanding  of  nuclear 
disturbed  environments  in  order  to  provide  reliable  predictions  of  the 
degradation  of  radar,  communications  and,  potentially,  advanced  weapons 
systems. 

The  experiment,  which  has  been  recently  upgraded  and  is  referred  to  as 

PHAROS  III,  illuminates  one  side  of  a  planar  target  with  a  beam  from  a 

neodymium  laser.  The  target,  usually  comprised  of  aluminum,  is  a  few  microns 

thick  and  the  laser  energy  is  typically  100  Joules.  The  laser  intensity  at 

13  2 

the  target  surface  is  about  10  Watts/cm  and  the  laser  wavelength  is  1 
micron.  During  the  laser  pulse  of  -  5  nanoseconds  FWHM,  the  front  side  of 

the  target  absorbs  the  laser  light  predominantly  by  inverse  bremsstrahlung 
and  heats  up  to  a  temperature  of  -  1  keV.  The  resulting  overpressure  leads 

to  a  forward  ablative  blovoff  with  velocities  of  -  10^  cm/sec  toward  the 
laser,  and  a  slower  backward  acceleration  of  the  rearside  of  the  target. 
When  a  background  gas,  typically  molecular  nitrogen,  is  present  in  the 
surrounding  chamber  at  pressures  ranging  from  0.5  to  5  Torr,  the  forward 
expanding  debris  forms  a  blast  wave  as  it  sweeps  up  the  background.  The 
blast  wave  reaches  1  cm  at  -  100  nanoseconds.  The  rearward  moving  debris 

also  forms  a  shock  wave  in  the  nitrogen,  but  its  precise  character  depends  on 
the  target  thickness.  This  pressure  regime  is  suitable  for  studying  strong 
to  moderate  collisional  coupling  between  the  debris  and  background  gas.  At 
lower  pressures,  a  kilogauss  magnetic  field  can  be  produced  around  the  target 
in  order  to  study  collisionless  coupling  mechanisms. 

As  part  of  the  theoretical  effort  to  understand  the  blast  wave  formation 

in  the  high  pressure  regime,  a  numerical  simulation  code  was  developed  within 

2 

the  Plasma  Physics  Division  .  This  simulation  code  modeled  the  dynamics  with 
a  Lagrangian  difference  scheme  in  one  spatial  dimension.  It  began  with  a 
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laser  pulse  impinging  on  a  cold,  aluminum,  planar  target.  The  target 
disassembly  occurred  by  heating  due  to  inverse  bremsstrahlung  absorption  of 
the  laser  beam.  While  the  target  remained  in  the  near  solid  phase,  the  ideal 
equation  of  state  was  corrected  to  include  the  zero  temperature  isotherm  of 
solid  aluminum.  The  subsequent  vaporization  and  hydrodynamic  expansion  of 
the  debris  and  swept-up  background  gas  vas  followed  out  to  ~  1  cm  assuming 
perfect  coupling,  i.e.,  a  contact  discontinuity  was  maintained  between  the 
debris  and  background.  The  nearly  radial  expansion  of  the  debris-air 
interface  was  accounted  for  by  employing  a  unique  oblate  spheroidal  geometry 
with  the  symmetry  axis  along  the  initial  laser  beam.  The  initial  target 
position  is  in  the  central  part  of  the  geometry  which  is  planar.  Separate 
ion  and  electron  temperatures  were  included  as  well  as  the  thermal  conduction 
from  each  species.  Magnetic  fields  were  not  included  for  these  high  pressure 
models.  In  order  to  follow  the  atomic  physics  starting  from  the  dense,  solid 
phase  to  the  low  density  gas  phase,  a  time-dependent  chemistry  model  was  used 
developed  including  all  the  ion  and  neutral  species  of  both  aluminum  and 
nitrogen.  Saha  ionization  equilibrium  conditions  are  maintained  while  the 
aluminum  is  near  solid  density.  After  vaporization,  both  the  debris  and 

background  chemistry  evolved  in  a  non-equilibrium  manner.  A  simple  radiation 
transport  scheme  was  developed  to  model  the  heating  of  the  backside  of  the 
target  due  to  x-rays  from  the  critical  surface  at  the  front  side. 

Some  of  the  conclusions  of  the  numerical  simulations  were  concerned  with 
the  nature  of  the  cavity  or  bubble  which  forms  behind  the  forward  moving 

blast  wave.  It  was  found  that  this  region  is  hotter  (T  -  oOeV)  and  more 

17-3  e 

rarefied  (n  -  2x10  cm  )  than  the  surrounding  blast  wave  (T  -  20eV  and  n 

18-3  e  e 

-  3x10  cm  ).  However,  the  cavity  is  highly  overionized  for  its 

temperature,  with  the  debris  species  predominantly  Al*  ,  Al**  ,  and  Al^  . 

The  results  for  the  blast  wave  region  are  in  good  agreement  with  the 

3 

experimental  results  ,  as  is  also  the  electron  density  in  the  cavity,  at 

4 

least  for  that  part  of  the  cavity  near  the  blast  wave  .  The  high  charge 
state  result  in  the  cavity  was  unexpected,  though  no  experimental 

determination  one  way  or  the  other  is  presently  known.  There  are  several 

simplifications  in  the  above  described  numerical  code  which  might  weaken  the 
conclusions.  First,  the  assumption  of  a  contact  discontinuity  between  the 
debris  and  air  is  not  strictly  valid.  •  •<’<=>n  nr  3  Tori  pressure.  The  aluminum 


and  nitrogen  will  mix  to  some  degree  and  subsequent  charge  exchange  between 
the  two  elements  could  lower  the  charge  state  of  aluminum  in  the  mixing 
region.  Second,  the  simplified  radiation  treatment  used  only  12  frequency 
bins,  did  not  specifically  account  for  resonant  line  emission,  and  the 
chemistry  model  did  not  include  self-absorption  within  each  grid  cell.  The 
neglect  of  these  radiative  effects  could  alter  the  cavity  temperature. 

In  order  to  improve  the  previous  calculations  and  solidify  the 
conclusions  we  have  upgraded  the  simulation  code  in  two  major  ways.  The  new 
code  is  referred  to  as  the  HANEX  code.  For  the  first  upgrade,  the  cross¬ 
streaming  of  the  aluminum  debris  through  the  nitrogen  background  is  followed 
approximately.  Since  the  basic  simulation  code  is  hydrodynamic,  it  models 
the  ion  and  electron  gases  as  a  Maxwellian  distribution  with  small  deviations 
due  to  temperature  gradients.  By  allowing  for  additional  deviations  to  the 
distribution  function  due  to  composition  and  partial  pressure  gradients, 
relative  drift  velocities  between  the  aluminum  ions,  nitrogen  ions,  and 
electrons  arise.  The  model  is  consistent  as  long  as  the  relative  drift 
velocities  are  small  compared  to  the  sound  speed.  The  resulting  mixing  by 
diffusion  of  the  highly  charged  aluminum  and  the  low  charged  nitrogen  in  the 
cavity  region  leads  to  charge  exchange  between  the  two  species.  The  theory 
behind  our  drift  calculation  and  charge  exchange  cross-sections  will  be 
reported  on  in  a  companion  report. 

The  second  major  upgrade  to  the  original  code  deals  with  radiation 
transport.  We  have  increased  the  frequency  resolution  of  the  radiation 
spectrum  from  12  to  166  frequency  bins.  This  permits  resolving  the  edges  of 
the  ground  states  of  all  aluminum  and  nitrogen  species  and  including  41 
resonance  transitions.  We  have  also  improved  the  calculation  of  the  emission 
and  absorption  processes.  The  theory  behind  our  second  upgrade,  along  with  a 
discussion  of  the  background  air  chemistry,  will  be  discussed  in  §11  of  this 
report.  The  processes  of  radiation  absorption,  subsequent  ionization,  and 
heating  in  the  background  plasma  are  addressed  ir  §111.  We  present  details 
of  these  processes  at  a  position  2  cm  away  from  the  target  and  over  the  time 
between  the  peak  of  the  laser  pulse  to  just  before  the  shock  front  sweeps 
over  this  position.  Also  model  spectra  from  the  simulation  runs  will  be 
presented  in  §111.  These  spectra  are  calculated  as  if  a  spectroscopic  camera 
viewed  the  experiment  perpendicular  to  the  laser  axis.  The  highly  ionized 
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emission  lines  from  resonant  transitions  which  appear  in  the  model  spectra 
are  indicative  of  the  highly  ionized  cavity.  Direct  comparison  with  observed 
spectra  in  the  UV  to  X-ray  energy  bandpass  should  provide  one  direct  test  of 
the  validity  of  our  simulation  models. 


II.  RADIATION  AND  PHOTO-CHEMISTRY 
A.  Emission,  Absorption,  and  Transport. 

The  treatment  of  radiation  in  the  HANEX  code  follows  the  same  basic 
transport  scheme  used  in  the  earlier  simulation  code;  the  radiation  intensity 
is  calculated  along  a  single  ray  according  to  the  fundamental  radiative 
transfer  equation, 


-  K  I 
V  \) 


Here  1^  is  the  radiation  intensity  in  ergs/cm  -sec-sd-Hz,  is  the 
absorption  coefficient  in  cm-1,  and  j  is  the  emission  coefficient  in 
ergs/cm  -sec-sd-Hz.  Each  of  these  quantities  are  carried  at  position  s  and 
frequency  v,  where  s  is  the  path  length  along  the  single  ray  parallel  to  the 
laser  axis  ar.d  passing  through  the  target  center. 

In  the  earlier  code  the  intensity,  absorption,  and  emission  coefficients 
were  averaged  over  12  frequency  bins.  In  the  HANEX  code,  however,  the 
resolution  of  the  radiation  spectrum  is  vastly  improved  by  employing  166 
frequencies  between  1  eV  and  10  keV.  The  radiation  transport  is  then 
calculated  frequency-by-frequency  and  the  absorbed  radiation  energy  is 
averaged  into  frequency  bins.  The  frequencies  ire  chosen  to  resolve  a  number 
of  spectral  features  composed  of  edges  and  lines.  Each  ground  state  of 
aluminum  and  nitrogen  leads  to  an  absorption  edge  for  the  outermost  valence 
electron  subshell,  giving  20  edges.  Each  edge  is  resolved  by  two 
frequencies,  one  slightly  above  and  the  other  slightly  below  the  edge. 
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Furthermore,  41  strong  resonance  lines  between  the  ground  and  excited  states 
are  included,  with  at  least  one  from  each  ionic  species.  The  specific  lines 
are  listed  by  ion  species  and  energy  in  Table  I.  (The  letters  in  parentheses 
refer  to  Fig.  7  discussed  in  section  Ill.b.)  Each  resonance  line  is  resolved 
by  three  frequencies,  one  at  line  center  and  the  other  two  in  the  line  wings. 
The  line  resolution  is  too  crude  for  determining  line  profiles  for  it  treats 
each  line  as  a  triangular  distribution,  but  the  resolution  is  sufficient  to 
account  for  the  energy  exchange  between  the  dynamics  and  radiation. 

The  plasma  emission  processes  include  free-free  bremsstrahlung,  free- 
bound  radiative  recombination,  and  bound-bound  emission  lines.  The 
absorption  processes  include  free-free,  bound-free,  bound-bound,  and  inner 
shell.  For  the  free -free  processes  the  absorption  coefficient  is 


Kv,ff  ‘  % 
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-hw/kT 
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where  e  is  the  elementary  charge,  mg  is  the  electron  mass,  c  is  the  speed  of 

light,  k  is  the  Boltzmann  constant,  T  is  the  electron  temperature,  n  is  the 

6  2  e 

electron  density,  n^  is  the  ion  density,  and  <Z  >  is  the  averaged  squared- 
charge  ,  i .  e . , 

n  7  2 

2  ni  Zi 

<Zz>  =  Z  -J — i- 
J  n. 

and  n^  is  the  number  density  of  ions  in  ionization  stage  j  with  charge  Zy 
From  Kirchoff's  law  the  corresponding  emission  coefficient  is 


jw,ff  =  3  J 
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The  free-free  Gaunt  factor  g was  fixed  at  unity  throughout  the  range  1  eV  < 
h\>  <  10  keV. 

For  the  bound-free  processes,  let  ^ ( j , k)  be  the  cross-section  for 
photo-absorption  from  ionization  stage  j  and  excited  level  k  to  the  next 
higher  ionization  stage  j+1.  Using  the  Milne  relation  to  connect  the  cross- 
section  for  photo-absorption  to  the  cross-section  for  induced  recombination 
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from  ion  stage  j+1  to  excited  level  k  of  ion  j,  the  photo-absorption  cross- 
section,  corrected  for  induced  recombination,  is  found  to  be 


where 
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Here  I.  ,  is  the  ionization  energy  of  excited  level  k  of  ionic  state  j,  the 
j  ,  x 

statistical  weight  of  the  same  level  is  g.  ,  ,  and  U.  is  the  partition 

J  *  *■  J 

function  for  the  ionic  state  j.  The  corresponding  emission  coefficient  due 
to  recombination  to  level  k  of  ion  state  j  is  given  by 


3  -hv/kT 

j  uf(J.k)  -  n.  .  a  (j,k)  0  e  e  .  (6) 

Jv,bf  J  j , k  v  J  2 

J  c 

Several  approximations  have  been  employed  in  these  equations  in  the 
HANEX  code.  First,  only  ground  states  are  considered  in  the  bound-free 
processes.  Second,  every  photo-absorption  cross-section,  a ,  is  assumed  to 
fall  off  as  l/v3  above  the  threshold  frequency.  Finally,  the  combination  of 
the  partition  functions  and  statistical  weights  in  0  is  taken  as  unity,  as  it 
would  be  in  local  thermodynamic  equilibrium  (LTE).  The  cross-sec t ion ,  a ^  ^f, 
at  threshold  is  the  maximum  of  two  values;  (i)  one  based  on  the  radiative 
recombination  coefficient  from  the  Milne  relation  (see  Appendix  1  of  Ref.  5), 
and  (ii)  the  other  based  on  the  inner  shell  opacity  data  discussed  in  the 
next  paragraph.  The  first  value  insures  detailed  balance  in  regions  of  high 
densi ty . 

In  addition  to  the  valence  shell  absorption  process  discussed  above  we 
included  the  inner  shell  opacity.  This  process  is  important  in  the  backside 
of  the  target  and  in  the  ambient  background  gas  where  the  flux  of  high  energy 
photons  formed  at  the  frontside  of  the  laser-heated  target  is  absorbed.  The 
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dependence  of  the  absorption  cross-section  on  energy  is  taken  from  Ref. 6  for 
the  neutral  species.  For  the  ionized  atoms  the  same  energy  dependence  is 
used  but  the  threshold  frequency  for  each  of  the  shells  shifts  to  higher 
energies.  These  shifts  are  given  by  Roothaan-Hartree-Fock  calculations  of 
the  atomic  level  binding  energies  in  Ref. 7. 

For  bound-bound  photo-absorption  from  the  ground  level  of  ion  stage  j  to 
excited  state  k,  the  cross-section  corrected  for  stimulated  emission,  is 

a  (j->k)  =  —  f..  *  (l  -  -J i kgj tiS  )  (7) 
v,bbVJ  mgc  jk  m  l  n^g.  J  K  ’ 

where  f.^  is  the  oscillator  strength  for  the  transition  and  <t>^  is  the  line 
profile  function.  Under  conditions  of  LTE  for  the  excited  states  the  term  in 
parentheses  reduces  to  l-expC-hu^/kT) .  The  corresponding  line  emission 
coefficient  is 
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J  v, bbv J  j  ,k  kj  4n  \ 


with  the  relation  between  the  Einstein  coefficients  being 
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Resolution  of  the  line  profile  would  require  many  frequencies  over  the 
Doppler  line  profile  Ve  have  simplified  the  analysis  by  adopting  a 
triangular  line  profile  which  vanishes  at  and  where  Am  is  the 
full  width  at  half  maximum  of  a  Doppler  profile.  At  line  center  the  profile 
function  is  then  <^(jk)  *  1/Am,  which  specifies  the  cross-section  for  the 
line  absorption  in  eqn.(7). 

The  cross  talk  between  our  time-dependent  chemistry  and  radiation 
transport  proceeds  in  the  following  manner.  After  a  hydro  timestep  6t  the 
chemistry  integrator  is  called.  The  chemistry  package  advances  the  species' 
abundances  over  St  and  simultaneously  calculates  for  each  grid  cell,  h,  the 
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frequency  and  time  integrated  emission  from  bremsstrahlung,  recombination, 
and  line  radiation.  To  complete  the  solution  for  the  timestep  St  in  our 
time-splitting  scheme,  the  radiation  package  is  called  last.  It  begins  by 
redistributing  the  integrated  emission  (from  the  chemistry  package)  into  a 
frequency  grid  for  each  cell  according  to  rhe  frequency  dependence  of  the 
different  j^'s  in  eqns.  (3),  (6),  and  (8).  Likewise  the  absorption 

coefficients  are  evaluated  from  eqns.  (2),  (4),  and  (7)  at  each  frequency, 

and  the  optical  depth  through  each  cell  h,  of  width  is  formed  as  x^  ^  1/2 
■  na^A^.  The  transport  equation  (1)  is  now  solved  for  the  radiation  energy 
at  the  cell  interfaces,  once  in  the  forward  direction  (  +  ) 
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and  once  in  the  rearward  direction  (-), 
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where  jv  },+  l/2  *s  t*ie  radiation  emitted  at  frequency  \>  from  all 

processes,  St  is  the  time  step,  Sv  is  the  frequency  interval  around  the 

frequency  \>,  and  *s  t*ie  v°lume  of  cell  h+1/2.  The  factor  1/2  in 

front  of  j  ,  -  means  that  one  half  of  the  photons  move  forward  and  the 

J v, h+1/2 

other  half  backward.  The  boundary  conditions  at  the  outermost  cell  interfaces 
are  that  no  radiation  energy  enters  the  system.  The  radiation  energy 
absorbed  in  cell  h  at  frequency  \>  is  given  by 
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The  last  term  in  this  equation  is  the  amount  of  energy  that  is  both  emitted 
and  then  re-absorbed  vithin  cell  The  chemistry  package  also  determines 

this  quantity,  but  under  slightly  different  conditions.  If  one  subtracts 
from  h+1/2  amount  of  radiated  energy,  AH,  which  the  chemistry  routine 
already  re-absorbed  vithin  cell  h,  the  remainder  gives  the  heating  of  the 

t 

electrons  in  cell  h;  H  =  H  .  .  ..  -  AH.,  ,  1/0.  The  amount  of 


The  amount  of 
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radiation  lost  from  the  system  is  given  by  the  sum  of  E+  at  the  forward  cell 

interface  and  E  at  the  rearvard  cell  interface.  If  molecular  nitrogen  is 

present  in  a  cell,  ^,+1/2  *s  as  discussed  below. 

Before  returning  to  the  next  hydro  and  chemistry  step  the  radiation 

routine  calculates  an  escape  probability  for  each  frequency  from 


1  _  e 'S>,h+l/2 
Pv’esc  =  \,h+l/2 

The  chemistry  routine  then  uses  this  probability  in  the  next  timestep  to 
determine  the  fraction  of  emitted  radiation  by  a  given  process  which  is  re¬ 
absorbed  within  the  same  cell  (  *  AH)  .  This  feature  takes  account  of  the 
re-absorption  of  radiated  energy  in  optically  thick  cells  and  therefore 
effectively  decreases  the  radiative  emission  leaving  a  cell. 

The  scheme  we  have  just  described  has  been  tested  for  consistency  by 
analyzing  a  high  density  lOeV  material.  The  resulting  spectrum  reproduced 
the  blackbody  limit  to  within  a  percent.  Furthermore,  the  simulation  runs 
conserved  total  radiative,  chemical,  internal  and  kinetic  energy  to  vithin 


The  above  scheme  does  have  the  limitation  that  it  cannot  handle  photo¬ 
pumping  and  photo-ionization  from  one  cell  to  another.  An  energy-consistent 
simulation  of  this  process  in  a  time-dependent,  non-equilibrium  chemistry, 
which  is  coupled  to  the  radiation  transport,  would  require  major  code 
development  and  testing.  At  present  those  photons  which  would  lead  to 
pumpimg  or  ionization  are  treated  the  same  as  the  remaining  photons,  i.e., 
they  are  deposited  as  electron  heating.  This  approximation  of  converting  all 
the  absorped  non-local  photons  to  pure  electron  heating  is  reasonable  in 
regions  near  local  thermodynamic  equilibrium  (LTE).  However,  in  the  cold, 
neutral  background  the  photo- ioniza r ion  due  to  the  non-local  photons 


predominantly  determines  the  state  of  the  gas.  The  special  treatment  of 
radiation  deposition  for  this  region  is  discussed  next. 


B.  Photo-Chemistry  of  the  Background  Gas 

The  phenomenology  in  the  region  ahead  of  the  shock  is  complex  and  far 
from  LTE.  In  this  section  we  describe  the  essential  approximations  adopted 
to  model  the  chemistry  for  this  region  in  the  HANEX  code.  During  the  laser 
pulse  the  ionization  and  heating  of  the  background  gas  is  driven  by  the 
absorption  of  the  UV  and  X-rays  from  the  laser  heated  target.  In  this 
deposition  phase  the  electron  temperature  in  the  background  gas  is  not  well 
defined  ,  i.e.,  the  electron  distribution  is  not  Maxwellian.  An  absorbed 
high  energy  photon  will  ionize  or  dissociatively  ionize  background 
molecules  and  create  energetic  electrons.  These  electrons  will  lose  most  of 
their  energy  by  further  ionization  and/or  electronic  excitation  of  molecules 
and  atoms.  Thermal  equilibration  of  the  fresh,  high  energy  electrons  with 
the  older,  lower  energy  electron  population  takes  place  primarily  through 
collisional  de-excitation  of  these  excited  states.  Radiation  losses  are 
negligible.  Further,  some  portion  of  the  electron  energy  will  excite 
vibrational  levels  of  the  N~  molecules.  But  the  relaxation  time  for 
vibrational  energy  is  long,  about  10  sec  for  a  5  Torr  background  pressure, 
so  vibrational  temperature  is  not  well  defined  either.  The  neutrals  and  ions 
will  be  heated  through  collisional  quenching  of  the  vibrationally  excited 
states.  Following  the  time  dependent  development  of  these  processes  in 
detail  is  beyond  the  scope  of  this  work.  In  order  to  approximate  the  state 
of  the  background  we  have  developed  a  simplified  model  to  account  for  these 
ionization,  dissociation  and  heating  processes. 

The  absorbed  photon  energy  is  apportioned  according  to  the  following 
prescription.  First,  the  cross  section  for  N?  is  assumed  to  be  twice  that  of 
N.  Then,  for  every  35  eV  of  the  absorbed  photon  energy  one  ion  is  created, 
either  N+  or  N2+,  depending  on  the  relative  abundance  of  N  and  N^.  In 
addition,  in  the  case  of  absorption  by  Nj  molecules,  for  every  55  ev  of 
absorbed  photon  energy  one  ^  is  dissociated.  The  remaining  energy  after 
ionization  or  dissociation  of  N.,  is  divided  evenly  between  the  heavy 
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particles  and  the  electrons.  For  photo-ionization  by  N  all  the  remaining 

energy  is  given  to  the  electrons.  Apportioning  some  of  the  energy  to  the 

heavy  particles  compensates  approximately  for  the  energy  fed  in  from 

vibrational  states,  which  is  not  tracked  in  detail. 

Time  dependent  chemistry  for  atomic  nitrogen  of  all  charge  states  has 

2 

been  described  in  a  previous  report  .  In  order  to  follow  the  behavior  of  the 
background  gas,  N2  chemistry  is  also  required.  The  following  two-body 
reactions  and  rate  coefficients  have  been  added  to  the  chemistry  model; 
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where  T  is  the  electron  temperature  in  eV  for  an  equivalent  Maxwellian 

e  2 
distribution.  In  the  last  reaction  the  products  are  really  N  and  N(  D). 

2 

However,  the  N(  D)  is  assumed  to  be  immediately  deexcited  by  electron 
collisions  with  the  2.38  eV/reaction  distributed  to  the  electron  gas. 

Although  this  prescription  may  seem  in  some  degree  arbitrary,  it  can  be 
justified  by  the  results  of  a  much  more  detailed  calculation  performed  by  two 

Q 

of  the  authors  using  a  similar  spectrum  for  the  deposited  radiation  .  The 
results  of  this  simple  model  for  fractional  ionization  and  dissociation,  as 
well  as  the  temperatures,  are  in  qualitative  agreement  with  the  detailed 
code. 


III.  DISCUSSION  OF  RESULTS 

We  begin  by  noting  that  the  total  radiation  loss  from  the  laser-target 
interaction  is  not  large.  Figure  la  displays  the  temporal  energy  budget  for 
our  standard  run,  vis. ,  a  100  J  laser  energy  shot  on  a  5.6  urn  thick  aluminum 
target  with  1.5  Torr  background  pressure  of  flv  Of  the  total  laser  energy, 
half  is  assumed  to  be  focused  and  absorbed  into  a  spot  o.  125  um  radius  on 
the  target.  For  a  laser  pulse  of  10  nsec  this  results  in  a  peak  intensity  of 
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-1.6x10  Watts/cm  .  (Note  that  the  125  um  is  the  radius  of  the  neck  in  the 
oblate  spheriodal  geometry  of  the  HANEX  code  pictured  in  Fig.l  of  Reference 
2.)  Figure  la  in  the  present  report  shows  that  the  time-integrated  energy 
lost  through  radiation  is  -20%  of  the  absorbed  laser  energy.  However, 
radiation  does  play  a  significant  role  in  transfering  energy  from  hotter 
regions  to  colder  ones.  Initially  the  laser  radiation  is  absorbed  by  inverse 
bremsstrahlung  in  a  thin  surface  layer  of  the  target.  Secondary  radiation 
from  this  hot,  dense  region  competes  with  electron  thermal  conduction  and  the 
impulse  shock  wave  in  transfering  energy  to  the  rear  side  of  the  target. 
Forward  propagating  radiation  from  this  same  hot  region  is  the  main  source  of 
energy  deposition  in  the  background  gas  during  the  laser  pulse.  At  later 
times,  radiation  from  the  frontside  blast  wave  becomes  important  both  for 
determining  the  condition  of  the  background  and  the  energy  loss  from  the 
system.  Figure  lb  displays  the  positions  of  the  frontside,  as  well  as  the 
rearside,  shockfronts  as  a  function  of  time.  In  the  next  subsection  we  will 
focus  our  attention  on  the  properties  of  the  background  gas  and  blast  wave. 
In  addition  to  energy  transfer,  radiation  losses  can  be  used  as  a  diagnostic 
tool.  In  the  following  subsection  we  will  discuss  how  UV  and  X-ray  spectra  of 
the  blowoff  plasma  may  be  used  to  discern  the  charge  species  in  the  cavity 
region  within  the  blast  wave. 

A.  Radiation  Deposition  in  the  Background  Gas. 

The  deposition  of  energy  in  the  gas  at  early  times  depends  on  the 
spectrum  of  the  radiation  from  the  target.  This  , however,  depends  on  the 
plasma  conditions  in  the  target.  In  Fig.  2  are  plots  of  the  electron  and  ion 
density,  velocity,  average  ion  charge,  and  electron  and  ion  temperatures  vs. 
radius  at  a  time  near  the  peak  of  the  laser  pulse  (5.48  nsec  after  the  pulse 
begins).  The  left  boundary  in  the  plots  coincides  with  the  rear  edge  of  the 
foil  and  the  extent  on  the  right  side  includes  only  part  of  the  blowoff 
region.  Three  distinct  regions  can  be  enumerated  in  these  plots.  A  cold 
dense  backside  with  some  rearward  expansion  near  the  foil  edge,  a  low 
density,  hot  blowoff,  and  a  transition  region  with  steep  gradients  in 
temperature,  density,  and  charge  state.  The  blowoff  region  is  composed  of 
high  charge  states  of  A1  and  will  emit  radiation  mostly  from  K-shell  lines 
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and  recombination.  Because  the  density  is  relatively  low,  especially  toward 
the  front,  most  of  the  generated  radiation  escapes  in  the  forward  direction. 
From  the  backside  through  the  transition  layer  the  Al  charge  states  increase, 
from  +3  to  +10.  The  radiation  produced  in  these  regions  will  therefore  be 
from  L-shell  lines,  recombination  and  bremsstrahlung.  Because  of  the  high 
densities  in  these  two  regions,  the  line  radiation  is  very  optically  thick 
and  almost  none  escapes.  Even  the  recombination  edges  have  a  optical  depth 
greater  than  one,  so  much  of  the  spectral  structure  will  be  absorbed.  The 
radiation  which  does  emerge  resembles  a  blackbody  at  a  temperature  near  that 
of  the  front  of  the  high  density  region.  This  radiation  will  pass  through 
the  blowoff  region  unattenuated  because  it  lies  below  the  K-absorption  edge. 
Thus  the  radiation  from  the  target  at  the  peak  of  the  laser  pulse  looks  like 
a  blackbody  near  40  eV  plus  a  high  energy  component  of  Al  K-shell  edges  and 
lines. 

A  plot  of  the  radiation  intensity  at  a  position  2  cm  from  the  target 
along  the  laser  axis,  for  the  case  described  above,  is  presented  in  Fig.  3a. 
Superimposed  on  the  target  spectrum,  one  can  clearly  see  the  absorption  edge 
due  to  the  intervening  N  and  Nj  between  the  target  and  the  above  position. 
The  energy  absorbed  by  the  background  at  this  2  cm  location  as  a  function  of 
the  photon  energy  is  shown  in  Fig.  3b.  The  spectral  regions  that  are  most 
important  for  absorption  at  this  time  of  the  peak  laser  pulse  lie  between  40 
and  200  eV  due  to  the  L-edge  of  nitrogen.  Another  absorption  peak  at  450  eV 
is  due  to  the  nitrogen  K-edge.  The  radiation  from  15  to  40  eV,  which 
actually  has  the  maximum  absorption  crossection,  has  already  been  absorbed  by 
atoms  nearer  the  target.  The  electrons  produced  by  photo-ionization  at  2  cm 
will  therefore  be  quite  energetic. 

After  the  peak  of  the  laser  pulse  the  plasma  rapidly  expands  and  the 
emitted  radiation  per  unit  volume  drops.  However,  the  amount  of  radiation 
that  reaches  2  cm  declines  more  slowly.  Although  the  emitting  target 
material  is  optically  thick  during  the  laser  pulse,  it  has  a  small  surface 
area.  After  the  material  expands  it  becomes  optically  thinner  and  the  volume 
emissivity  specifies  the  radiation.  The  large  volume  of  the  expanded  plasma 
maintains  the  amount  of  radiation  energy  incident  at  2  cm.  As  the  blowoff 
region  expands  into  the  nitrogen  background,  line  radiation  from  the  blast 
wave  becomes  an  important  source  of  energy  deposition  in  the  more  distant 
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background.  In  Figs.  4a  and  4b  the  incident  and  absorbed  spectrum, 
respectively,  are  shown  at  a  time  60  nsec  after  the  peak  of  the  laser  pulse. 
At  this  time  both  spectra  are  dominated  by  line  radiation  from  the  blast  wave 
as  is  evidenced  by  the  series  of  sharp  lines  near  30  -  60  eV.  The  fraction 
of  the  total  energy  incident  at  2  cm  in  each  of  four  spectral  intervals  is 
plotted  as  a  function  of  time  in  Fig.  5a.  The  shift  from  energetic  photons 
at  early  times  to  less  energetic  ones  at  late  time  can  be  seen  in  this 
figure.  In  Fig.  5b  the  fraction  of  the  total  energy  absorbed  at  2  cm  is 
shown  for  the  same  intervals.  The  most  energetic  interval  contains  photons 
above  667  eV.  This  radiation  is  from  the  K-shell  of  aluminum  and  originates 
in  the  blowoff  region.  The  radiation  originating  mostly  from  the  backside 
region  is  divided  into  two  spectral  intervals,  153  eV  to  398  eV,  and  398  eV 
to  667  eV.  The  lowest  interval,  photon  energy  less  than  153  eV,  is 
predominately  from  the  cold  target  regions  at  early  times  and  from  line 
radiation  from  the  blast  wave  at  later  times. 

The  effect  of  the  deposition  of  the  radiation  on  the  background  gas 
depends  on  the  energy  absorbed  per  atom.  In  Fig.  6a  the  time  integrated  and 
the  rate  of  energy  absorbed  per  atom  at  2  cm  is  shown.  The  nearly  constant 
rate  of  energy  absorption  after  20  nsec  is  due  to  the  fact  that  the  decrease 
in  the  radiative  output  from  the  shock  front  per  unit  area  is  compensated  by 
a  decrease  in  the  distance  to  the  shock  as  it  moves  outward  and  closer  to  2 
cm.  The  densities  of  ^  N,  N2+,  N+  and  the  electrons  ng  are  presented  vs. 
time  in  Fig.  6b.  We  note  that  the  fraction  of  ionization  is  only  a  few 
percent  during  the  laser  plulse  but  rises  to  nearly  15%  just  before  the  time 
that  the  shock  overruns  the  position  at  2  cm. 
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B.  Model  Spectra  of  the  Blowoff  Plasma. 

Inclusion  of  the  physics  of  radiation  emission  and  absorption  in  a 
simulation  code  is  important  for  modeling  the  dynamical  evolution  of  the 
laser-target  interaction  and  blowoff.  In  addition,  a  numerical  method  for 
radiation  transport  can  create  model  spectra  which  can  be  compared  with 
experiment.  It  often  happens  that  spec troscopic  observations  provide  the 
most  detailed  data  of  a  plasma  experiment. 
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To  produce  model  spectra  we  have  employed  the  previously  described 
radiation  transport  scheme  with  several  modifications  in  a  postproccesser . 
The  first  modification  is  a  change  in  the  transport  direction.  In  the  hydro¬ 
chemistry-radiation  code  the  radiation  was  transported  along  a  ray  parallel 
to  the  laser  beam  axis.  This  approach  successfully  accounts  for  the  transfer 
of  radiation  energy  from  an  inner  plasma  region  to  an  outer  one.  In  the 
laser  experiment  all  spectroscopic  observations  are  performed  with  the  camera 
line-of-sight  perpendicular  to  the  laser  axis.  To  compare  with  these 
observations  the  postprocessor  solved  the  radiation  transport  equation  (1) 
along  a  set  of  rays  perpendicular  to  the  laser  axis.  Each  ray,  along  which 
the  transport  equation  is  to  be  solved,  begins  at  the  far  side  of  the  plasma 
emitting  region.  After  passing  across  the  laser  axis  and  through  the  near 
side  of  the  emitting  plasma,  the  ray  encounters  an  extended  region  (10  cm 
long)  of  neutral  The  addition  of  this  absorbing  region  surrounding  the 
blowoff  plasma  is  the  second  modification.  Transport  through  it  models  the 
fact  that  between  a  spectroscopic  camera  and  the  experiment  lies  an  extended 
region  of  nitrogen  at  the  background  pressure.  Finally,  we  assume  in  the 
postprocessor  that  the  plasma  blowoff  is  spherical  about  the  laser  target. 
This  modification  from  the  geometry  of  the  hydrodynamic  code,  which  is  oblate 
spheroidal,  permits  a  great  simplification  in  computing  the  intensity  along 
line-of-sight  rays  which  pass  above  or  below  the  laser  axis.  Viewing  the 
blowoff  plasma  perpendicular  to  the  laser  axis,  the  emergent  intensity  along 
all  rays  which  pass  at  a  distance  r  from  the  target  is  equal  to  the  solution 
of  the  transport  equation  for  the  ray  which  crosses  the  laser  axis  at 
position  r.  Geometrical  considerations  then  allow  the  calculation  of  an 
approximate  model  spectra  at  any  spatial  position. 

As  an  example,  we  study  the  spectra  arising  from  a  laser  produced  plasma 
for  a  5  Torr  background  pressure  and  20  J  laser  energy  on  a  4.6  micron  thick 
aluminum  target  of  radius  125  microns.  Artificial  masks  were  used  so  that 
only  the  forward  moving  blowoff  plasma  would  contribute  to  the  model  spectra. 
Specifically,  if  the  target  and  blowoff  plasma  are  enclosed  by  a  large 
circle,  the  back  side  of  the  circle  away  from  the  laser  beam  was  blocked  off. 
Further,  a  disc  of  radius  0.25  cm  was  centered  on  the  initial  position  of  the 
target  to  block  out  emission  from  the  dense  target  material.  To  maximize  the 
emission  intensity,  the  model  spectra  was  time  integrated  over  the  first  100 


nsec  of  the  run.  Figure  7a  shows  the  resulting  model  spectrum.  The  low 

energy  continuum  is  due  to  bremsstrahlung  and  the  high  energy  jagged 

continuum  results  from  recombination  edges  of  nitrogen  and  aluminum.  The 

large  gap  in  the  continuum  between  15  eV  and  -50  eV  is  due  to  the  inner  shell 

absorption  by  the  10  cm  long  path  of  extended  nitrogen.  The  strong  emission 

lines  are  labeled  by  letters  according  to  Table  I.  The  model  spectrum 

includes  emission  from  the  blast  wave  formed  in  the  background  nitrogen  gas 

as  well  as  emission  from  the  cavity  within  the  blast  wave.  This  cavity  is 

composed  of  aluminum  and  nitrogen.  The  latter  has  been  mixed  in  by 

diffusic  The  high  energy  lines,  particularly  those  above  -100  eV,  are 

emitted  from  the  cavity.  The  parent  ionic  species  of  these  lines  from  both 

aluminum  and  nitrogen  reveal  a  highly  charged  region.  The  lower  energy  lines 

arise  from  the  blast  front.  In  the  simulation  for  this  run  the  electron 

temperature  is  -20  eV  in  the  blast  wave,  but  rises  to  -70  eV  in  the  cavity. 

Hence,  as  has  been  noted  in  the  Introduction,  the  cavity  is  overionized  for 

its  temperature.  This  situation  is  a  remnant  of  the  early-time  freezing  out 

of  the  charge  state;  the  initially  hot  target  material  cools  by  expansion  at 

a  faster  rate  than  the  recombination  rate.  Early  in  time  the  blast  wave  was 

also  much  hotter  and  the  swept-up  nitrogen  was  thereby  ionized  up.  Due  to 

the  difficulties  in  directly  measuring  the  density  and  temperature  of  the 

frontside  cavity  region,  we  suggest  that  experimental  spectra  in  the  UV  and 

X-ray  domain  should  be  quite  revealing  in  light  of  the  predicted  spectra  of 

Fig.  7a.  For  comparison,  a  similar  model  spectra  of  the  rear  side  is 

presented  in  Fig.  7b.  In  this  case  the  masks  were  turned  over  to  block  out 

the  frontside  and  target  emission.  At  late  times  the  cavity  stretches  from 

the  frontside  blastwave,  through  the  initial  target  position,  and  backwards 

to  near  the  rearward  moving  shock.  Hence  some  of  the  moderately  high 

ionization  stages  appear  both  in  the  front  and  backside  spectra;  for 

instance,  lines  i, j ,m,n,q,r, t,and  u.  Other  lines  from  moderately  high 

ionization  stages  of  aluminum,  such  as  h,k,o,and  p,  appear  on  the  backside 

but  not  the  frontside  spectrum.  The  striking  difference  between  the  backside 

and  frontside  spectra  is  the  strong  lines  from  low  ionization  stages  of 

aluminum  (lines  a  and  b)  on  the  rearside  spectrum  only,  in  contrast  with  the 

•*"11.  •*"  1  ^ 

cluster  of  lines  (v-z)  from  Al  and  Al  on  the  frontside  spectrum  only. 
These  features  from  the  spectra  reflect  the  details  of  the  simulation  run: 
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the  frontside  cavity  is  more  highly  ionized  than  the  rcarside  one,  while  the 
rearward  moving  blast  wave  is  denser  and  cooler  than  its  frontside 
counterpart . 


IV.  SUMMARY  AND  DISCUSSION 

In  this  report  we  have  presented  a  review  of  the  multi-frequency 
radiation  transport  scheme  used  in  the  NRL  HANEX  simulation  code.  Likewise 
we  have  discussed  the  model  incorporated  in  HANEX  for  the  absorption, 
ionization,  excitation,  and  heating  of  the  background  nitrogen  resulting  from 
the  deposition  of  the  high  energy  target  and  blast  wave  radiation.  Results 
were  presented  for  the  evolution  of  the  background  nitrogen  from  the  time  of 
the  peak  in  the  laser  pulse  to  60  nsec  later.  Model  spectra  of  the  cavity 
and  blast  wave  region  indicate  strong  lines  from  highly  ionized  aluminum  and 
nitrogen.  This  result  suggests  that  UV  and  X-ray  spectroscopy  of  the  plasma 
blowoff  could  be  used  to  investigate  the  nature  of  the  cavity  and  verify  our 
long-standing  prediction  of  a  hot  (50-70  eV) ,  overionized  (<Z>~5-10)  cavity 
within  the  blast  wave. 

The  transport  scheme  and  deposition  model  described  in  section  II 
employed  several  simplifying  assumptions.  First,  the  boundaries  of  the 
oblate  spheroidal  geometry,  so  useful  for  the  dynamic  expansion,  impose  the 
artificial  condition  that  the  radiation  is  constrained  to  diverge  through  a 
cone  of  45  degrees  half-angle.  At  late  times  when  the  radiation  is 
predominately  from  the  blast  wave  this  condition  will  produce  a  small  error 
since  any  deposition  occurs  nearby.  But  at  early  times  this  condition  will 
over  estimate  the  flux  at  a  distant  background  position.  Fortunately,  the 
time  integrated  radiation  escaping  the  target  at  early  times  is  low  compared 
to  the  final  value  of  about  15  to  20%.  Second,  the  treatment  of  line 
transport  is  very  approximate  and  leads  to  an  overestimate  of  the  line 
absorption  when  the  lines  are  moderately  thick  as  in  the  dense  target 
regions;  no  allowance  is  made  for  line  radiation  to  escape  through  broad 
wings  of  the  line  profile.  Third,  the  assumption  that  all  the  energy 
absorbed  in  a  cell  is  put  into  heating  of  the  electrons  will  not  cause  errors 
if  the  cell  is  near  equilibrium,  as  in  the  dense  target  region  and  the  blast 
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wave.  However,  when  the  target  material  falls  to  low  density  (the  cavity) 
the  plasma  is  far  from  equilibrium.  This  material  does  not  absorb  much 
energy  because  the  absorption  cross-section  is  near  zero  for  the  low  energy 
photons  that  dominate  our  spectrum.  In  this  region  a  more  careful  treatment 
of  absorbed  radiation  would  somewhat  decrease  the  electron  temperature.  The 
problem  of  depositing  radiation  in  the  cold  neutral  background  gas  has  only 
been  approximated.  The  energy  absorbed  in  the  nitrogen  K-shell  will  lead  to 
further,  though  lover  energy,  radiation  and  therefore  not  be  deposited 
locally.  This  process  could  effect  the  temperatures  of  neutrals  and  ions,  as 
well  as  the  electrons.  It  is  believed  that  the  prescription  used  for  the 
deposition  process  should  provide  values  for  densities  to  within  a  factor  of 
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TABLE  I 

TRANSITION  ENERGY  (in  eV)  OF  TRANSPORTED  RESONANCE  LINES 


A1  I 

4.01 

5.23 

5.78 

A1  II 

7.42(b) 

A1  III 

6.67(a) 

17.81 

Al  IV 

77.44(o) 

94. 95(p) 

A1  V 

44. 34(k) 

Al  VI 

40.03(h) 

Al  VII 

34.92 

Al  VIII 

37.95 

49.74(1) 

Al  IX 

43.73(j) 

31.8 

41 . 0( i ) 

Al  X 

37.24(g) 

238. 5(q) 

Al  XI 

22.30 

256. 6(r) 

Al  XII 

1598. 0(v) 

1869. 0(x) 

1964. 0(y) 

Al  XIII 

1724. 0(v) 

2042. 0(z) 

N  I 

10.33 

N  II 

11.43(d) 

13.53(f) 

19.22 

N  III 

12.51(e) 

16.23 

18.08 

N  IV 

16.20 

50. 15(m) 

N  V 

10.00(c) 

59.24(n) 

N  VI 

430.67(s) 

497. 94( t) 

N  VII 

500. 0(u) 

23.22 

33.12 
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Fig.  7(a)-A  synthetic  spectrum  of  the  forward  moving  plasma  blowoff  and  blast 
wave  integrated  in  time  (up  to  100  nsec)  and  in  space.  The  rearside  and 
target  region  where  blocked  off.  Prominent  lines  are  marked  by  letters 
corresponding  to  those  transitions  listed  in  Table  I.  The  initial 
conditions  are  for  a  20  J  laser  shot  with  a  5  Torr  background  pressure. 
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Fig.  7(b)-A  synthetic  spectrum  for  the  rearward  moving  plasma.  In  this  case 
the  frontside  and  target  have  been  blocked  off. 
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MS  D4  3  8 

LOS  ALAMOS,  NM  87545 
0 1 C Y  ATTN  S.P.  GARY 
01 C Y  ATTN  J.  BOROVSKY 

SANDIA  LABORATORIES 
P.O.  BOX  5800 
ALBUQUERQUE,  NM  871 1 5 

0 1  C  Y  ATTN  DOC  CON  FOR  W.  BROWN 
0 1 C Y  ATTN  DOC  CON  FOR  A. 

THORNBROUGH 

01CY  ATTN  DOC  CON  FOR  T.  WRIGHT 

0 1 C Y  ATTN  DOC  CON  FOR  D.  DAHLGREN 

0 1 C Y  ATTN  DOC  CON  FOR  3141 

01 CY  ATTN  DOC  CON  FOR  SPACE  PROJECT 
DIV 

SANDIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.O.  BOX  969 
LIVERMORE,  CA  94550 

01 CY  ATTN  DOC  CON  FOR  B.  MURPHEY 

0 1 C  Y  ATTN  DOC  CON  FOR  T.  COOK 


OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  DC  20545 

0 1 C Y  ATTN  DOC  CON  DR.  YO  SONG 

NATIONAL  OCEANIC  4  ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 
0 1 C Y  ATTN  R.  GRUBB 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORPORATION 

P.O.  BOX  92957 

LOS  ANGELES,  CA  90009 


01  CY 

ATTN 

I. 

GARFUNKEL 

0  1  C  Y 

ATTN 

T. 

SALMI 

01  CY 

ATTN 

V  . 

JOSEPHSON 

0  1  C  Y 

ATTN 

S. 

BOWER 

01  CY 

ATTN 

D. 

OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CCRP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

0 1 C  Y  ATTN  RADIO  SCIENCES 

AUSTIN  RESEARCH  ASSOC.,  INC. 

1 901  RUTLAND  DRIVE 
AUSTIN,  TX  78758 

0 1  C  Y  ATTN  L.  SLOAN 
01 CY  ATTN  R.  THOMPSON 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.O.  BOX  983 
BERKELEY,  CA  94701 

0 1 C  Y  ATTN  J.  WORKMAN 
0 1 C  Y  ATTN  C.  P RETT  IE 
0 1 C  Y  ATTN  S.  BRECHT 

BOEING  COMPANY,  THE 
P.O.  BOX  3707 
SEATTLE,  WA  98124 

0  1  C  Y  ATTN  G.  KEISTER 
0 1  C  Y  ATTN  D.  MURRAY 
01CY  ATTN  G.  HALL 
0 1 C  Y  ATTN  J.  KENNEY 

CHARLES  STARK  DRAPER  LABORATORY,  INC 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 
0 1 C  Y  ATTN  D.3i  COX 
0 1 C  Y  ATTN  J.  P.  GILMORE 
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COMSAT  LABORATORIES 
22300  COMSAT  DRIVE 
CLARKSBURG,  MD  20871 
0  1  C  Y  ATTN  G.  HYDE 


CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  1  *4  850 

0 1 C  Y  ATTN  D.T.  FARLEY,  JR. 


ELECTROSPACE  SYSTEMS,  INC. 
BOX  1359 

RICHARDSON,  TX  75080 
01CY  ATTN  H.  LOGSTON 
01 CY  ATTN  SECURITY  (PAUL 


PHILLIPS) 


EOS  TECHNOLOGIES,  INC. 

606  Wilshire  Blvd. 

Santa  Monica,  CA  90*401 
01  C  Y  ATTN  C.B.  GABBARD 
0 1 C  Y  ATTN  R.  LELEVIER 


GTE  SYLVANIA,  INC  . 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM  ,  MA  021  9*4 

0  1  C  Y  ATTN  DICK  STEINHOF 


HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD ,  MA  01730 
0  1  C  Y  ATTN  DONALD 


HANSEN 


OF 


ILLINOIS,  UNIVERSITY 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 
0 1 C  Y  ATTN  K.  YEH 


INSTITUTE  FOR  DEFENSE  ANALYSES 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA ,  VA  22311 
0 1 C  Y  ATTN  J.M.  AEIN 
0 1 C  Y  ATTN  ERNEST  BAUER 
0 1 C  Y  ATTN  HANS  WOLFARD 
0 1 C  Y  ATTN  JOEL  BENGSTON 


I  NT  L  TEL  &  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY ,  NJ  07110 

0 1 C  Y  ATTN  TECHNICAL  LIBRARY 


JAYCOR 

11011  TORREYANA  ROAD 

P:0; -BOX  8515*4 

SAN  DIEGO ,  CA  92138 

0 1 C  Y  ATTN  J.L:  SPERLING 


GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

(ALL  CLASS  ATTN:'  SECURITY  OFFICER) 

0 1 C  Y  ATTN  T.N.  DAVIS  (UNCLASS  ONLY) 
0 1 C  Y  ATTN  NEAL  BROWN  (UNCLASS  ONLY) 


JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

0 1 C  Y  ATTN  DOCUMENT  LIBRARIAN 
0 1 C Y  ATTN  THOMAS  POTEMRA 
01 C  Y  ATTN  JOHN  DASSOULAS 


KAMAN  SCIENCES  CORP 
P.0.  BOX  7*463 

COLORADO  SPRINGS,  CO  80933 
0  1  C  Y  ATTN  T.  MEAGHER 


KAMAN  TEMPO-CENTER  FOR  ADVANCED 
STUDIES 

816  STATE  STREET  (P.0  DRAWER  QQ ) 
SANTA  BARBARA,  CA  931  02 
0 1 C  Y  ATTN  DASIAC  ' 

0 1 C  Y  ATTN  WARREN  S.  KNAPP 
0 1 C Y  ATTN  WILLIAM  MCNAMARA 
0  1  C  Y  ATTN  B.  GAMBILL 


LINKABIT  CORP 
1  0*4  53  ROSELLE 
SAN  DIEGO ,  CA  92121 

0 1 C  Y  ATTN  IRWIN' JACOBS 


LOCKHEED  MISSILES  i  SPACE  CO, 
P.O.  BOX  50*4 
SUNNYVALE,  CA  9*4  088 

0 1 C  Y  ATTN  DEPT  60-12 
0 1 C  Y  ATTN  D.R.  CHURCHILL 
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LOCKHEED  MISSILES  &  SPACE 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 


ATTN  MARTIN  WALT  DEPT  52-12 
ATTN  W.L.  IMHOF  DEPT  52-12 
ATTN  RICHARD  G.  JOHNSON 
DEPT  52-12 

ATTN  J.B.  CLADIS  DEPT  52-12 


MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  V A  22101 

0  1  C  Y  ATTN  W.  HALL 
0 1 C  Y  ATTN  W.  FOSTER 


MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 

0 1 C Y  ATTN  R.  HEFFNER 


MCDONNEL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 


ATTN  N.  HARRIS 
ATTN  J.  MOULE 
ATTN  GEORGE  MROZ 
ATTN  W.  OLSON 
ATTN  R.W.  HALPRIN 
ATTN  TECHNICAL 
LIBRARY  SERVICES 


MISSION  RESEARCH  CORPOR 
735  STATE  STREET 
SANTA  BARBARA,  CA  931  01 


CORPORATION 


MISSION  RESEARCH  CORP. 
1720  RANDOLPH  ROAD,  S.E. 
ALBUQUERQUE,  NM  87106 
0 1 C  Y  R.  STELLINGWERF 
0 1  C  Y  M.  ALME 
01 C  Y  L.  WRIGHT 


PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

01 C  Y  ATTN  E.C.  FIELD,  JR. 


PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

01 C  Y  ATTN  IONOSPHERIC  RESEARCH 


PHOTOMETRICS,  INC. 

4  ARROW  DRIVE 
WOBURN,  MA  01801 

01  CY  ATTN'  IRVING  L  . 


KOFSKY 


PHYSICAL  DYNAMICS,  INC. 
P.O.  BOX  3027 
BELLEVUE,  WA  98009 

0 1 C Y  ATTN  E.J.  FREMOUW 


0  1 C  Y 

ATTN 

P.  FISCHER 

PHYSICAL  DYNAMICS, 

INC  . 

01  CY 

ATTN 

W.F.  CREVIER 

P.O.  BOX  10367 

01CY 

ATTN 

STEVEN  L.  GUTSCHE 

OAKLAND,  CA  94610 

01  C  Y 

ATTN 

R.  BOGUSCH 

ATTN  A.  THOMSON 

0  1  C  Y 

ATTN 

R.  HENDRICK 

01  CY 

ATTN 

RALPH  KILB 

R  4  D  ASSOCIATES 

0  1  C  Y 

ATTN 

DAVE  SOWLE 

P.O.  BOX  9695 

01CY 

ATTN 

F.  FAJEN 

MARINA  DEL  REY ,  CA 

90291 

0 1  C  Y 

ATTN 

M.  SCHEIBE 

0 1 C  Y  ATTN  FORREST  GILMORE 

01  CY 

ATTN 

CONRAD  L.  LONGMIRE 

0 1 C  Y  ATTN  WILLIAM  3.  WRIGHT 

0  1  C  Y 

ATTN 

B.  WHITE 

0 1 C  Y  ATTN  WILLIAM  J.  KARZAS 

0 1  C  Y 

ATTN 

R.  STAGAT 

0 1 C  Y  ATTN  H.  ORY 

0  ICY 

ATTN 

D.  KNEPP 

0 1 C  Y  ATTN  C.  MACDONALD 

01  CY 

ATTN 

C.  RINO 

0  1  C  Y  ATTN  BRIAN 

LAMB 

MORGAN  GROVER 


RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 
■  01CY  ATTN  BARBARA  ADAMS 
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RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  U2nd  STREET 
NEW  YORK,  NY  10036 

0  1  C Y  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS 

INTERNATIONAL  INCORPORATED 
1 1 50  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 

0 1 C Y  ATTN  LEWIS  M.  LINSON 
0 1 C Y  ATTN  DANIEL  A.  HAMLIN 
01 CY  ATTN  E.  FRIEMAN 
0  1  C Y  ATTN  E . A .  STRAKER 
01CY  ATTN  CURTIS  A.  SMITH 

SCIENCE  APPLICATIONS 

INTERNATIONAL  CORPORATION 
1 71 0  GOODRIDGE  DR. 

MCLEAN,  VA  22102 
0  1  C  Y  J.  COCKAYNE 
01  CY  E.  HYMAN 

SRI  INTERNATIONAL 


TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD  ,  MA  01730 

0  1  C Y  ATTN  W.P.  80 QU 1ST 


TRW  DEFENSE  4  SPACE 
ONE  SPACE  PARK 
REDONDO  BEACH, 

ATTN  R. 

ATTN  S. 

ATTN  D. 

ATTN  D/ 


\v\  wy wy-»  vwrji 

£ 


CRCUP 


0  1  C  Y 
01  CY 
0  1  C  Y 
Cl  C  Y 


CA  90273 
K.  PLE3UCH 
ALTSCHULER 
DEE 

STOCKWELL 


SNTF/1 575 


V ISIDYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON ,  MA  01803 
0 1 C  Y  ATTN  W.  REIDY 
0  1  C  Y  ATTN  J.  CARPENTER 
0 1 C  Y  ATTN  C.  HUMPHREY 

UNIVERSITY  OF  PITTSBURGH 
P ITTSBURGH ,  PA  15213 

01  CY'  ATTN:  N.  ZABUSKY 


din;: 
u. 


333  RAVENSWOOD  AVENUE 

ANNrM'Oi 

MENLO  PARK,  CA  9«025 

01  CY 

ATTN 

J.  CASPER 

01  CY 

ATTN 

DONALD  NEILSON 

01  CY 

ATTN 

ALAN  BURNS 

Code  2628 

0 1  C  Y 

ATTN 

G.  SMITH 

01  CY 

ATTN 

R.  TSUNODA 

0  1  C  Y 

ATTN 

DAVID  A.  JOHNSON 

0 1  C  Y 

ATTN 

WALTER  G.  CHESNUT 

Code  1220 

0 1 C  Y 

ATTN 

CHARLES  L.  RINO 

01  CY 

ATTN 

WALTER  JAYE 

0  1  C  Y 

ATTN 

J.  VICKREY 

Records 

01  CY 

ATTN 

RAY  L.  LEADABRAND 

0  1  C  Y 

ATTN 

G.  CARPENTER 

01  CY 

ATTN 

0.  PRICE 

0  1  C  Y 

ATTN 

R.  LIVINGSTON 

01  CY 

ATTN 

V.  GONZALES 

0  1  C  Y 

ATTN 

D.  MCDANIEL 

~j'r  AR  C 

2  cop ie: 


22  copies 


1  copy 


1  CODV 
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